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This article illustrates how very small deviations from the Maxwellian expo- 
nential tail, while leaving unchanged bulk quantities, can yield dramatic effects 
on fusion reaction rates and discuss several mechanisms that can cause such 
deviations. 

Fusion reactions are the fundamental energy source of stars and play im- 
portant roles in most astrophysical contexts. Since the beginning of quantum 
mechanics, basic questions were addressed such as how nuclear reactions occur 
in stellar plasmas at temperatures of few keV (1 keV « 11.6 x 10^ °K) against 
Coulomb barriers of several MeV and what reactions or reaction networks dom- 
inate the energy production. It was soon realized that detailed answers to such 
questions involved not only good measurements or quantum mechanical un- 
derstanding of the relevant fusion cross sections, but also the use of statistical 
physics for describing energy and momentum distributions of the ions and their 
screening 

Gamow understood that reacting nuclei penetrate Coulomb barriers by means 
of the quantum tunnel effect and Bethe successfully proposed the CNO and then 
the pp cycle as candidates for the stellar energy production: this description has 
been directly confirmed by several terrestrial experiments that have detected 
neutrinos produced by pp and CNO reactions in the solar core ^5]. 

In the past only few authors (e.g., d'E. Atkinson, Kacharov, Clayton, Haubold) 
examined critically the energy distribution and proposed that such distribution 
could deviate from the Maxwellian form. In fact, it is commonly accepted that 
main-sequence stars like the Sun have a core, i.e. an electron nuclear plasma, 
where the ion velocity distribution is Maxwellian. In the following, we first 
discuss why even tiny deviations from the Maxwellian distribution can have 
important consequences and then what can originate such deviations. 



Thermonuclear reaction in plasmas and distribu- 
tion tails 



In a gas with ni (712) particles of type 1 (2) per cubic centimeter and relative 
velocity v, the reaction rate r (the number of reactions per unit volume and 
unit time) is given by 

r = {1 + Si2)~'^ nin2{av} , (1) 

where a = <^{v) is the nuclear cross section of the reaction. The reaction rate 
per particle pair is defined as the thermal average 

f{v)avdv, (2) 







where the particle distribution function f{v) is a local function of the tempera- 
ture PI- 

Therefore, the reaction rate per particle pair (av) is determined by the spe- 
cific cross section and by the velocity distribution function of the reacting parti- 
cles. When no energy barrier is present and far from resonances, cross sections 
do not depend strongly on the energy. Most of the contribution to {av) comes 
from particles with energy of the order of kT, and the dependence on the spe- 
cific form of f{v) is weak. The same is true for bulk properties that receive 
comparable contributions from all particles: e.g., the equation of state. 

The situation is very different in the presence of a Coulomb barrier, when 
the reacting particles are charged, as in the fusion reactions that power stars [7j- 
The penetration of large Coulomb barriers {Za/r is of the order of thousands 
in units of kT when r is a typical nuclear radius) is a classically forbidden 
quantum effect. The penetration probability is proportional to the Gamow 
factor exp [— Eq /E] , where the Gamow energy Eq = 2/ic^(ZiZ2a7r)^, a is 
the fine structure constant, /i is the reduced mass, and are the charges 
of the ions. The cross section is exponentially small for E ^ Eq and grows 
extremely fast with the energy; therefore, one usually defines the astrophysical 
S factor, whose energy dependence is weaker 

a(E)^^e-V^ . (3) 

The two factors in the integrand in Eq. (0) that carry most of the energy de- 
pendence are the Maxwellian distribution cx e'~'^^^^^\ which is exponentially 
suppressed for E ^ fcT, and the penetration factor e~ v^^Ab^ which is expo- 
nentially suppressed for E <C Eq. Contributions to the rate come only from an 
intermediate region (Gamow peak) around the temperature-dependent energy 
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which is called the most effective energy, since most of the reacting particles 
have energies close to Eq. 

Figure Ogives a pictorial demonstration of how the Gamow peak originates 
and how different reactions select different parts of the distribution tail and can 
be used to probe it. 
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In the upper panel (a) the exponentiaUy decreasing function (thin black 
curve) is the Maxwellian factor; the rapidly growing function (dash blue curve) 
is the penetration factor (for graphical reason multiplied by 10^) of one of the 
most important reactions in the Sun, ^He + ^He ^Re + 2p [Eg — 11.83 MeV), 
which corresponds to a most effective energy Eq — (EcikT)^ /Ay^^ — 17.036 kT 
for kT = 1.293 keV = 11.6 x lO*' °K; the product of the two functions (Gamow 
peak) is the thick red curve. Note that the Gamow peak, and therefore the 
rate, is very small (it has been multiplied by an additional 10® factor to make it 
visible on the same scale of the other curves), since at the most effective energy 
Eq both the cross section and the number of particles are exponentially small. 
At this point is important to remark that the area under the Maxwellian curve 
for energies within the Gamow peak (the energy window indicated by the red 
band) is of the order of 0.1% of the total area: only a few particles in the tail 
of the distributions contribute to the fusion rate. 

The fact that the penetration factor effectively selects particles in the tail 
of the distribution is the more dramatic the larger the charge of the reacting 
ions: for the p+ ^^N — > ^^O + 7 (the leading reaction of the CNO cycle, which 
dominates the energy production in main-sequence stars larger or older than 
the Sun) the contributing particles are few in a million. 

The effect on the Gamow peak when increasing the charges of the reacting 
nuclei is shown in the lower panel (b) of Figure Q The green, red, and blue 
curves show the Gamow peak multiplied by 10^, 10^^, and 10^^, respectively, 
for three fundamental reactions in main-sequence stars: p + p d + v + e'^ 
{Eg = 493 keV), ^He + ^He ^ ^Kc + 2p {Eg = 11.83 MeV), and p + 

^^O + 7 {Eg = 45.09 MeV)). It is immediately evident that the larger the 
charges of the ions the higher is the energy of the particles that contribute to 
the rate, the (much) lower is the peak and, therefore, the (much) smaller is the 
rate. In fact the maximum of the Gamow peak Eq cx E^^ cx (^1^2)^/"^ and its 
hight is proportional to exp {—3Eo/kT). 

A convenient parametrization of deviations from the Maxwell distribution is 
the deformed g-exponential: 



which naturally appears in Tsallis' formulation of statistical mechanics i4 . This 
particular deformation of the exponential has the advantage of describing both 
longer tails (for q > 1) and cut-off tails (for q < 1), while reproducing the 
exponential in the limit g — > 1. 

Figure 13 shows the effect of substituting to the Maxwellian distribution 
exp{—E/kT) the distribution Nqexp^{—E/kT), where Nq is the normalization 
factor that conserves the total number of particles. We show the effect for the 
p + reaction and the values (a) q = 1 ± 0.002 and {h) q = 1 ± 0.01. This 
reaction determines the rate of energy production from the CNO cycle, domi- 
nant at older stages and, therefore, also determines the time of the exit from 
the main sequence. Black curves refer to the exponential {q = 1), red curves 
refer to the cut-off (g < 1) exponential, and green curves refer to the longer-tail 
{q > 1) exponential. 

Note that all exponentials have been multiplied times a huge 10^ factor for 
emphasizing their tiny differences: these values of q produce functions that are 
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almost indistinguishable from the exponential unless one looks very far in the 
tail. 

One can make several remarks: 

• Gamow peaks are shifted towards higher energies, when the distribution 
has a tail longer than the exponential (green curves, q > 1); they are 
shifted towards lower energies for cut-off exponentials (red curves, q < 1); 

• the effect is larger the larger the deviation from the exponential (the larger 

• the peaks (and the rates) become correspondingly higher for 5 > 1 and 
smaller for 5 < 1; 

• the effect on the rate is already large for |1 — = 0.002, it becomes huge 
(more than a factor of 10) for \l — q\ = 0.01; note that the green (red) peak 
in lower panel (b) of figure[21(|l — q\ — 0.01) has been divided (multiplied) 
by five to make them fit on the same scale! 

These deviations should be carefully estimated, since reliable calculations 
of nuclear reaction rates in stellar interiors is fundamental for a quantitative 
understanding of the structure and evolution of stars. In fact, while the overall 
stellar structure is rather robust, changes of some of the rates even by few 
percent can produce detectable discrepancies, when precise measurements are 
possible, e.g., in the case of the solar photon and neutrino luminosity, and 
mechanical eigenfrequencies In quasi-stellar objects like Jupiter deviations 
could be even larger and explain their excess energy . 

As already shown in the recent past, very small deviations from Maxwellian 
momentum distribution do not modify the properties of stellar core and are in 
agreement with the helioscismology constraints , but may affect the evaluation 
of the nuclear fusion rates that may be enhanced or depleted, depending on 
superdiffusion or subdiffusion property of the particles |Zj . 

Deviations from Maxwellian distribution 

Normal stellar matter, such as the one in the Sun, is non-degenerate, i.e., quan- 
tum effects are small (in fact, they are small for electrons and completely neg- 
ligible for ions), it is non-relativistic, and it is in good thermodynamical equi- 
librium. On this ground, the particle velocity distribution is almost universally 
taken to be a Maxwell-Boltzmann (MB) distribution. 

Concerning the thermodynamical equilibrium, main sequence stars are more 
precisely in a stationary state where the luminosity equals the heat production 
rate. This metastable state has a long life-time, of the order of the star life- 
time, and it ends when the nuclear fuel is burned out. In addition, the quasi- 
equilibrium is only local, since the temperature decreases from core to surface. 
However, nuclear reactions are often, but not always, sufhciently slow on the 
scale of thermal and mechanical exchanges and take place on such a small scale 
that spatial and temporal deviations from equilibrium can be neglected to a 
very good first approximation. 

At least in one limit the MB distribution can be rigorously derived: systems 
that are dilute in the appropriate variables, whose residual interaction is small 
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compared to the one-body energies. In spite of the fact that the effects of the 
residual interaction cannot be neglected (the electron screening factor is a well- 
known example of correction due to the astrophysical plasma environment) at 
zero order the many-body correlations can be neglected and the stellar interior 
can be studied in this dilute limit. In this limit the velocity distribution is the 
Maxwellian one. 

However, one should keep in mind that derivations of the ubiquitous Maxwell- 
Boltzmann distribution are based on several assumptions 7 . In a kinetic ap- 
proach, one assumes (1) that the collision time be much smaller than the mean 
time between collisions, (2) that the interaction be sufficiently local, (3) that 
the velocities of two particles at the same point are not correlated (Boltzmann's 
Stosszahlansatz), and (4) that energy is locally conserved when using only the 
degrees of freedom of the colliding particles (no significant amount of energy 
is transferred to collective variables and fields). In the equilibrium statistical 
mechanics approach, one uses the assumption that the velocity probabilities of 
different particles are independent, corresponding to (3), and that the total en- 
ergy of the system could be expressed as the sum of a term quadratic in the 
momentum of the particle and independent of the other variables, and a term 
independent of momentum, but if (1) and (2) are not valid the resulting ef- 
fective two-body interaction is non-local and depends on the momentum and 
energy of the particles. Finally, even when the one-particle energy distribution 
is Maxwellian, additional assumptions about correlations between particles are 
necessary to deduce that the relative- velocity distribution, which is the relevant 
quantity for rate calculations, is also Maxwellian. 

In the following we give arguments and mechanisms that lead to distribution 
functions that are different from the MB one in a stationary state. 

Correlations between particles, so that the probability distribution of the 
system is not described by the product of independent probabilities of the com- 
ponents, are in general responsible for such more general distributions. The 
specific microscopic mechanisms that generate these correlations depend on the 
particular system and there exist many approaches to derive the relevant dis- 
tributions. 

In an approach that uses the Fokker-Planck equation, which takes into ac- 
count the average effects of the environment through the drift J{p) and diffusion 
D{p) coefficients, stationary solutions different from the Maxwell distribution 
(e.g. Druyvenstein or Tsallis like distributions) can be obtained, when J{p) 
and D{p) include powers of p higher than the lowest order The presence of 
higher powers of p, i.e., higher derivative terms, can be interpreted as a signal of 
non-locality in the Fokker-Planck equation. We stress that these distributions 
are stationary (stable or metastable) and what counts to decide the distribu- 
tion is the type of collisions between ions and the dependence on momentum 
of the elastic collisional cross sections (Coulomb, screened Coulomb, enforced 
Coulomb, among others), or the presence of ion-ion correlations [HI- 

The presence of random fields (e.g., distributions of random electric micro- 
fields or, in general, of random forces) introduces in the kinetic equations factors 
whose effect is to enhance or to deplete the high-momentum tail of the distri- 
bution function 7 . 

Because of the many-body nature of the effective forces, which makes the 
collisions not independent, the distributions of the relevant degrees of freedom 
observed, e.g., the ones selected by a fusion reaction, can be different from the 
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distributions of the quasi-particles that describe the plasma. In addition the 
plasma makes effective interactions time dependent (memory effects) and non- 
local. These effects depend strongly on the energy of the selected particles and 
on the coUisional frequency. 

One important and clear example of this last point is given by the fact that 
many processes, such as nuclear fusion itself, depend on momentum rather than 
on energy. This distinction is important because, due to plasma many-body 
effects, an uncertainty relation holds between momentum and energy l1 Oj . Even 
when the energy distribution maintains its Maxwellian expression, the momen- 
tum distribution can be different in the high energy tail. In fact, this quantum 
uncertainty effect (not Heisenberg uncertainty) between energy £ and momen- 
tum p, caused by the many-body collisions and described by the Kadanoff-Baym 
equation, implies an energy-momentum distribution of the form 

/q(£,p) = -n{E)5^{£,p) (6) 

TT 

with 

where l^^{£,p) is the mass operator of the one-particle Green function. After 
integrating in d£ the product of fq^S^p) and the Maxwellian energy distribu- 
tion, we obtain a momentum distribution with an enhanced high-momentum 
tail. Although this approach produces a deviation from MB distribution, the 
state represented by /g (p) is an equilibrium state m . The Maxwellian 
distribution is recovered in the limit when Sj{£,p) becomes a S function with a 
sharp correspondence between momentum and energy. 

Distributions different from the Maxwellian one can also be obtained ax- 
iomatically from non-standard, but mathematically consistent, versions of sta- 
tistical mechanics that use entropies different from the Boltzmann-Gibbs one ^ 

US]. 

We have argued that it is not sufficient to know that the Maxwellian dis- 
tribution is a very good approximation to the particle distribution. We must 
be sure that there are no corrections to a very high accuracy, when studying 
reactions that are highly sensitive to the tail of the distribution, such as fusion 
reactions between charged ions. Several mechanisms have been outlined (others 
need to be studied) that can produce small, but important deviations in the tail 
of the distribution. 
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Figure 1: The Gamow peak and energy selection. In the upper panel (a) 
the exponential thin black curve is the Maxwellian distribution. The rapidly 
increasing dash blue curve shows the behavior of the penetration factor N x 
exp {-^JEg/E) for the solar reaction ^He + ^He ''He + 2p {Eg = 11.83 MeV 
and N — 10^). The red thick curve shows the product of the two curves (Gamow 
peak) times 10*. The horizontal red band indicates the energy range of the 
reacting particles. The lower panel (b) shows how different reactions select 
different windows of particle energies. The Gamow peak (energy window) moves 
to higher energies going from p + p ^ d + i^ + e+ {Eg — 493 keV, green), to 
^He + ^He ^ ''He + 2p {Eg = 11.83 MeV, red), and p+^^^ ^ ^^O + 7 
{Eg — 45.09 MeV, blue). Correspondingly, the peaks become (much) lower; 
note that the three curves have been multiplied times 10^, lO'^^, and 10^^, 
respectively, to make them visible on the same scale. 
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Figure 2: Effects of tiny changes in the tail of the distribution for the reaction 
p + ^^N. Using the parametrization of Eq. © , the upper panel (a) shows the 
effect of taking q — 1.002 (green) and q — 0.998 (red), while the lower panel (b) 
shows the effect of g = 1.01 (green) and q — 0.99 (red); black curves correspond 
to the usual Maxwell distribution {q — 1). In the lower panel (b) the green (red) 
Gamow peak (thick lines) is divided (multiplied) by an additional factor of five. 
Note that thin curves (exponentials and q-exponentials) have been multiplied 
by 10^ to enounce their tiny differences. 
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